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A new paradigm emerging in the description of magnetic-sensitive radical-ion-pair recombination 
reactions is presented. Tliis paradigm is founded on the realization that the recombination process 
of radical-ion pairs is a continuous quantum measurement. In the regime of low magnetic fields 
we describe the appearance of purely quantum phenomena, that fade away as the magnetic field 
increases. We will analyze the magnetic sensitivity of these reactions under this perspective and 
bridge the new full quantum theory with the existing classical reaction theory applicable at high 
magnetic fields. Based on the quantum theory of recombination we will then explain experimental 
observations incompatible with classical reaction theory, in particular the effect of deuteration on 
the magnetic sensitivity of radical-ion pair recombination yields. 



I. INTRODUCTION 

Radical-ion pairs have been at the focus of chemical, 
biological and physical research for the last few decades, 
mostly because of their relevance to photosynthesis pj] . It 
is a chain of electron-transfer reactions that converts 
the photon energy absorbed by pigment molecules in the 
photosynthetic reaction center to energy capable of driv- 
ing further chemical reactions of biological importance. 
The radical-ion pairs formed in this intricate photochem- 
ical machinery further complicate the picture, since the 
spin interactions of their unpaired electrons make the 
whole process magnetic-sensitive. The importance of the 
magnetic sensitivity of photosynthetic reactions [1, 0, Q 
was realized early on and intensely studied theoretically 
as well as experimentally. The timing of photosynthetic 
reactions apparently rests on the complicated interplay 
between non-radiative electron-transfer rates, spin-state 
mixing rates and spin-selective recombination rates, to 
name a few of the participating processes. It is therefore 
not surprising that theories and experimental techniques 
from all of the aforementioned disciplines have been sum- 
moned to address this problem. 

Parallel to these developments, it was early on real- 
ized @, i?]] that the magnetic sensitivity of radical- ion pair 
(RIP) recombination reactions might be at the basis of 
the biochemical magnetic compass used by several species 
for navigation in earth's magnetic field [8| . After the de- 
velopment of extensive theoretical models M. [lol [Til [l3 | 
and several supporting experiments [H, [l4 Il5l|. the so- 
called RIP magnetoreception mechanism is now estab- 
lished as forming the basis of the biochemical compass, 
at least of migratory avian species. 

The magnetic sensitivity of RIP reactions 1(| can be 
easily understood: the RIP is formed from a photoexcited 
donor-acceptor dyad molecule (DA) in the singlet spin 
state. Magnetic interactions with the external magnetic 
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field and hyperfine interactions with the molecule's nu- 
clear spins induce the so-called singlet-triplet (S-T) mix- 
ing. Since the recombination of the RIP is spin-selective, 
i.e. the singlet ^(D+A~) RIP can recombine to the neu- 
tral DA molecule, whereas the triplet "^(D+A^) recom- 
bines, if at all, to a triplet metastable excited state """DA, 
it is seen that the magnetic field can directly influence the 
end-result of the photoexcitation-recombination cycle. 

What has gone unnoticed so far, however, is the fact 
that the spin-selective charge recombination process is 
fundamentally a quantum measurement |17j . The tun- 
neling of the RIP into a vibrational excited state of the 
DA molecule is a scattering process the outcome of which 
depends on the RIP 'a spin state. In other words, the in- 
termolecular recombination dynamics continuously inter- 
rogate the radical-ion pair's spin state. This interroga- 
tion can be readily described by quantum measurement 
theory [T^ and leads to purely quantum phenomena such 
as the quantum Zeno effect [1^ . The theoretical descrip- 
tion of magnetic-sensitive RIP recombination processes 
was so far based on phenomenological classical reaction 
theory [2^ which masked the quantum effects that are 
dominant in the regime of low-magnetic fields. This is 
the reason that classical theories were rather successful, 
i.e. most measurements were performed at high mag- 
netic fields, where classical reaction theories are valid. 
Several puzzling experimental observations at low fields 
were however left unexplained. As will be shown in the 
following, the quantum Zeno effect crucially affects the 
relevant time constants of the RIP recombination pro- 
cess, and naturally leads to an understanding of several 
experimental findings that could not be reconciled with 
classical RIP reaction theories. In this article we will 
present the detailed quantum theory of the RIP recom- 
bination process and elucidate the fundamentally differ- 
ent interpretation of the low-magnetic-field effects that 
emerges. We will outline the connection of the full quan- 
tum theory to be presented with the classical reaction 
theories. We will also provide an explanation of experi- 
mental findings that have challenged the classical theory, 
in particular measurements regarding the effect of deuter- 
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FIG. 1: Quantum measurement scheme of the radical-ion 
pair charge recombination process. The coherent evolution of 
the RIP induced by the magnetic Hamiltonian is interrupted 
by the measurements performed by the charge-recombination 
process. At some random time the measurement is completed 
and a definite measurement result is available, so that the re- 
combination can proceed accordingly. 



ation on the magnetic-sensitivity of RIP recombination 
yields. 



II. THE QUANTUM DESCRIPTION OF RIP 
RECOMBINATION 

In Figure 1 we diagrammatically describe the evolu- 
tion of the RIP recombination. The RIP is at time t — Q 
created in the singlet state ^(D+A^), which coherently 
mixes with the triplet ''"(D'^A^) under the unitary ac- 
tion of the magnetic hamiltonian Ti.. At the same time, 
the charge recombination is performing a quantum mea- 
surement of the singlet- and triplet-state projection op- 
erators Qs and Qt, with the measurement rate being ks 
and kx, respectively. At some point the measurement 
is over. If the measurement result is qs = ^ (obviously 
meaning that qr = 0), the RIP can tunnel to an excited 
vibrational state of the neutral DA molecule, ^(DA)*, 
and from there decay to the DA ground state at a rate 
75. In the event that the measurement produces the re- 
sult qs — (obviously meaning that qr = 1), the RIP 
similarly ends up at a metastable triplet state of the DA 
molecule, """DA, and this happens at a rate 7t. In this 
work we will not deal with the decay rates 75 and 7t, 
since they not are relevant for the magnetic dynamics. 

Quantum measurement theory [l^ describes the effect 
of a continuous measurement of a system observable q 



performed at a rate k on the system's time evolution by 
the Liouville equation: 



^ = - k[q, [q,p\] 



(1) 



where p is the quantum system's density matrix and the 
first term describes the unitary evolution due to a Hamil- 
tonian interaction Ti.. In the case of the RIP recombina- 
tion process, the measured system observables are Qs 
and Qt- Therefore, the RIP's density matrix evolution 
equation will be 



dp 
'dt 



-i[n,p] - ks[Qs, [Qs,p]] - kriQT, [Qt,p]] (2) 



Measurements performed on a coherently evolving quan- 
tum system always cause decoherence. That is, ^ de- 
scribes the evolution of an initially pure state of the RIP 
to a mixed state, during which the coherence induced 
by the Hamiltonian H is also dissipated. This is due to 
the measurement-induced back-action. The RIP mag- 
netic Hamiltonian to be considered is comprised of the 
Zeeman interaction of the two unpaired electrons with 
the external magnetic field and by the hyperfine inter- 
actions with the molecule's nuclear magnetic moments. 
We will consider for simplicity only two nuclei, with nu- 
clear spins Ii and I2, and different hyperfine couplings 
(in frequency units) Ai and A2 to each of the electrons. 
We will consider isotropic hyperfine interactions, as we 
are not interested in describing directional effects [llj in 
this work. The magnetic Hamiltonian thus reads 



n = hi + h2 



(3) 



where 



(4) 



is the interaction Hamiltonian for each electron and j 
1,2. The Larmor frequency is 



(5) 



and throughout we are using units such that h = 1. For 
earth's field of ~ 0.5 G, w = 0.7 ^s~^. In the following, 
u; will be the unit of frequency, i.e. we will normalize 
all rates by tu. The normalization is performed by uj 
and not by the hyperfine frequency scale ai (or 02) for 
the following reason. If the mechanism is to serve for 
the measurement of the magnetic field uj, with precision 
6lu, then from Heisenberg's time-energy uncertainty rela- 
tion follows that in a measurement time t, the achievable 
precision Slu is such that tSuj ~ 1. The precision Slo is 
understood to be the width of the measured distribution 
of frequencies lo. Since the measurement time r is de- 
termined by the decay rate A of the RIP S-T mixing, 
T ^ 1/A, the relative width of the measured distribu- 
tion is Slo/u ~ X/lo. Therefore Heisenberg's uncertainty 
poses the requirement that X/uj < 1, if the mechanism is 
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to be sensitive to the magnetic field itself. If however, the 
signal-to-noise of the actual detection process is s, than 
the measurement sensitivity is enhanced to X/suj. For 
example, if the magnetic-field change is detected by the 
change in the recombined molecules yield and the latter 
can be measured at the level of 1%, then Suj/uj ^ 0.01 if 
A/w ~ 1. It is obvious that the condition X/uj < 1 can 
be more easily satisfied at high fields, since the classical 
recombination rates k (and hence the decay rates A) are 
usually much larger than 1 /is~^. In the following, we 
are going to show how the quantum Zeno effect allows 
the existence of long S-T lifetimes, so that the condition 
A/w 1 is satisfied even if » 1 ^s~^. 

We will now derive some basic properties of ([2|). We 
denote hy S ^ Tr{pQs} and T = Tr{pQT} the probabil- 
ity to find the RIP in the singlet and triplet spin state, 
respectively. Due to the fact that Qs and Qt are orthog- 
onal projection operators, we have = Qs, = Qt, 
QsQt = and Qs + Qt = 1- Therefore, using the 
identity Tr{[A, B]} = 0, it follows from ^ that 

^--»Tr{[Qs,H]p} (6) 

^ = -zTr{[QT,H]p} (7) 

and 

Tr{p} = S + T = 1 (8) 

Does the fact that Trjp} = 1 mean that the RIP never 
recombines? Not at all. Every time the RIP's spin state 
approaches the point S — 1 (or T = 1), there is a non- 
zero probability of singlet- or triplet channel charge re- 
combination. If the RIP recombines, it ceases to exist 
in the Hilbert space defined by the eigenstates of 7i, and 
together with is ceases the meaning of the density matrix 
p and its normalization. This is formally described by a 
quantum stochastic process [H, taking place in an 
enlarged Hilbert space. A brief analysis has appeared in 
[l3| and a more detailed presentation of the RIP's spin 
state time evolution will be presented in a forthcoming 
manuscript. Such a description requires the complete 
solution of the density matric equation ([2]). In the fol- 
lowing we will present some general arguments that do 
not necessitate such a solution. 



III. EIGENVALUE SPECTRUM OF THE 
LIOUVILLE EQUATION 

Since the nuclear spin multiplicity of the problem we 
are considering is (2/i + l)(2/2 + 1), and the electron 
spin multiplicity is 4, the density matrix is a square A^- 
dimcnsional matrix, where A = 4(2/i-t-l)(2/2+l). Equa- 
tion thus represents a set of A^^ coupled first-order 
differential equations. We can rewrite them in the form 
dR/dt — AR, where R — (pn pi2 ... PnnY' ^ column 
vector of dimension N"^ containing all density matrix el- 
ements and A is a A^^-dimensional square matrix. The 
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FIG. 2: The decay rates and mixing frequencies of ([2]) are 
sfiown as a function of ks/io, for fcr = Q.2ks, Ii ~ I2 = 1/2, 
and a2 ~ 2ai — Sui. 

eigenvalues of A are playing a central role in the consid- 
erations of this work, since they are governing the time 
evolution of the RIP's spin state. These eigenvalues are 
found by numerically diagonalizing the matrix A, and in 
general are complex numbers of the form —A + iil, where 
A > is termed the decay rate and O the mixing fre- 
quency of this particular eigen-mode. There are n = N"^ 
such eigenmodes. For example, if /i = 1/2 and /2 = 1, 
then n — 576. It is easily seen that n grows exponen- 
tially with the number of participating spins, and the 
numerical diagonalization quickly becomes a formidable 
problem. This is the reason we limit the calculations to 
the case of just two nuclear spins. 

A result of such a diagonalization is shown in Figure 
2 for typical values of the problem's parameters. We as- 
sume that kx < ks, so that the recombination rate scale 
is determined by ks- We will distinguish two physically 
distinct regimes: (i) ks/uj ^ 1, and (ii) ks/uj ^ 1. We 
term the first regime classical, since it is in this case that 
classical reaction theories are valid, whereas the second is 
termed quantum, because it is in this parameter regime 
that the quantum Zeno effect is at play. Before exam- 
ining the magnetic sensitivity in these two regimes, we 
note the following. Since the time evolution of the den- 
sity matrix is described by n eigenmodes, every density 
matrix element can be written as 

n 

p.,{i)=Y.^e^''^^^^^\ (9) 
1=1 

where the coefficients cf} are determined from the initial 
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conditions. Therefore 

n 

y 1=1 

where pij = {i\p\j) and \i) with i — 1, ...,n are the basis 
states, and 

withqji — (jIQsI*}- The reahty of ((5s(t)) is based on the 
appearance of positive and negative mixing frequencies, 
as shown in Figure 2B. 

A. Classical Regime 

We here study the magnetic sensitivity of RIP reac- 
tions when feg/cj < 1, i.e. in the regime of small recombi- 
nation rates and/or large magnetic fields. This is termed 
the classical regime. We first note that in this regime 
the decay rates A scale proportionally to ks, A ^ ks- 
Secondly, the mixing frequencies are constant, i.e. inde- 
pendent of the recombination rate ks- Hence all n terms 
contribute to {Qsit)), i.e. risiow = n in (H^. In this 
case, the magnetic sensitivity of {Qs{t)) stems from the 
dependence of qji on uj, i.e. a change in lo will result in 
a change of Ai. As can be seen in Figure 3B, the matrix 
elements qij and hence Ai are varying strongly for fields 
uj/a> I until the very high field regime w/a ^ 1, where 
the magnetic sensitivity drops to zero. In this high-field 
regime, the T± branches of the energy eigenvalues have 
separated away and there is small mixing only between 
To and S, as can be seen in Figure 3 A. The above con- 
siderations form the classical interpretation of magnetic 
field effects on RIP recombination reactions, that have 
already been analyzed in detail [§, (and references 
therein) . 

B. Quantum Regime 

The situation changes dramatically as the ratio ks/^ 
increases. Indeed, it is seen from Figure 3C that for low 
magnetic fields, to /a <^ 1, a 10% change in uj produces 
a mere 2% change in qij. According to classical theory, 
this is further suppressed by the ratio w/fcg, resulting in a 
negligible change of {Qs{t)). However, significant recom- 
bination yield changes are indeed observed experimen- 
tally. The reason is that in this regime quantum effects 
become dominant, fundamentally altering the nature of 
the magnetic sensitivity of the recombination process. As 
seen in Figure 2A, it is in this regime that the decay rates 
split into two branches, the fast (upper branch) and slow 
(lower branch) decay rates. The latter are a manifesta- 
tion of the quantum Zeno effect fl^ , the physical signif- 
icance of has been already described elsewhere [1^ . We 
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FIG. 3: (A) Energy eigenvalues of the magnetic Hamiltonian 
for nuclear spins /i = J2 = 1/2 and hyperfine couplings ai = 
a and 02 = 2a. (B) The corresponding matrix elements qij. 
In this case the Qs-matrix is 16-dimensionaI, so there are 
256 matrix elements qij of Qs in the basis formed by the 
eigenvectors of 7i. It has been checked that the eigenvalues of 
the Qs-matrix in this representation are and 1. (C) Zoom 
in at low magnetic fields, in order to better display the change 
of Qij in this regime. 

will here elaborate on the new interpretation of magnetic 
sensitivity of RIP reactions in this quantum regime. Let 
us suppose that some of the decay rates A/ are slow, i.e. 
A; <C ri/, and some other are fast, i.e. A; ^ ili- Let us 
further suppose that there are rigiow of the former and 
hence rifast = n — risiow of the latter. Thus nfast terms in 
PU)) will quickly decay away, leaving the risiow terms to 
determine the long-time behavior of {Qs{t))- 

"alow 

(Qs(0)*»i/A.„. = E e(-^'+'^^')*^, (12) 
1=1 

It thus follows that in this regime it is mainly the change 
of the number rigiow with the magnetic field lu that is 
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FIG. 4: (A) Decay rates only the lower branch is shown) as 
a function of u for ks ~ 15, kr = 0.2fcs, nuclear spins I\ — 
I2 = 1/2 and hyperfine couplings ai = 1.5 and 02 = 3. The 
dashed red line is the X — u line. The magnetic sensitivity 
derives from the changing number of decay rates crossing this 
line. 



responsible for the magnetic sensitivity of {Qs{t)). This 
number is a strongly changing function of uj as shown in 
Figure 4. 



IV. PHENOMENOLOGICAL LIOUVILLE 
EQUATIONS 



The phenomenological density matrix evolution equa- 
tions used so far are 

^ = -i[n, p] - ks{QsP + pQs) - kT{QTP + pQt) (13) 



The decay rates of the eigenvalue spectrum of p3)) are 
shown in Figure 5. It is seen that all the decay rates 
scale in proportion to ks, i.e. the quantum Zeno effect 
has disappeared. This is expected, because equation 
by design forces the normalization of the density matrix, 
Tr{p}, to an exponential decay: 



dt 



= ^ks{Qs)-kT{QT) 



(14) 



There is an analogous situation in atomic physics exper- 
iments: transit time broadening [23]. When an atomic 
coherence (for example a ground state spin coherence) 
of an atomic vapor is probed by a laser, the atoms stay 
in the interaction volume (defined by the laser beam) 
for a limited time Ttr, determined by the atomic veloc- 
ity distribution. The measured spin coherence lifetime 
will thus be limited by Tt,., which will add to the spin- 
resonance width by l/rtr- This is the transit-time broad- 
ening, which would be described by an equation similar 

to iir 



dTr{p} Tr{p} 



FIG. 5: Decay rates of classical reaction theory Liouville 
equation (|13p . The calculation was done for kr ~ 0.2fcs, 
Ii — I2 = 1/2, Hi = 1.5aj and a2 — 3uj. 



This equation phenomenologically describes the disap- 
pearance of atoms with a time constant Ttr. The situ- 
ation with radical-ion pairs, however, is fundamentally 
different, since in this case there is a selective quantum 
measurement continuously going on. An attempt to arti- 
ficially introduce the crucial term QspQs existing in Q 
and missing from (jl3p was made in [25)] based on phe- 
nomenological considerations, without any further con- 
sequences on the theoretical description of RIP recom- 
bination dynamics. All theoretical treatments, summa- 
rized in [2fl|], have thereafter used the phenomenological 
density matrix equation (|13|) . 



dt 



Ttr 



(15) 



V. EXPLANATION OF EXPERIMENTAL 
OBSERVATIONS INCONSISTENT WITH 
CLASSICAL REACTION THEORY 



We will now apply the interpretation of the low-field 
magnetic sensitivity of RIP recombination reaction previ- 
ously put forward to address experimental observations 
that seem paradoxical within the framework of classi- 
cal reaction theory. Unusually long time constants that 
have been observed in RIP reactions several times have 
been explained elsewhere ll^- Recently observed [iH] 
magnetic sensitivity of RIP reaction rates which is in 
stark contrast with classical reaction theory has been 
shown 27] to be fully consistent with the quantum the- 
ory. Several experimental observations relevant to avian 
magnetoreception based on magnetic-sensitive RIP reac- 
tions have also been accounted for [26|. Here we focus 
on some other puzzling observations, namely the effect of 
the rip's deuteration on the magnetic sensitivity of the 
RIP recombination. 
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A. Magnetic Sensitivity and Deuteration 

The effects of deuteration on the magnetic sensitivity 
of RIP recombination reactions have early on [2^ sig- 
nalled the inadequacies of classical reaction theories, as 
has been noted in J^. As stated in [2^, "The identi- 
cal triplet quantum yields in H and D samples does not 
support the generally accepted idea that hyperfine inter- 
actions are responsible for spin rephasing in the PSTO"*"!" 
radical pair" . Indeed, it is rather straightforward to de- 
rive, based on the classical spin-state mixing description, 
the result of a change of hyperfine couplings. The triplet 
recombination yield is crudely expected '291 to scale as 
a/ kg, where a is the dominant hyperfine coupling and fcg 
the singlet recombination rate. Any change in a should 
thus lead to a proportional change in the yield, contrary 
to the 1979 observations. 

It is rather straightforward to explain this apparent 
contradiction by use of the eigenvalue spectrum of the 
quantum evolution equation ([2]). Indeed, if we calculate 
the eigenvalue spectrum for a molecule with h ^ I2 ^ 
1/2 and hyperfine couplings ai = 1.5uj and 02 = 2.8ai, 
we find that in the quantum regime {ks /to 3> 1) the num- 
ber of decay rates A such that X/ui < 1 are 58.2% of the 
total (256). If we repeat the calculation for Ii = 1/2, 
I2 — 1, ai — 1.5lu and 02 — 0.86ai, i.e. we replaced 
the second nucleus, which before was assumed to be a 
hydrogen, with a deuterium, we find that the previous 
percentage changes to 62.5% (the total number of eigen- 
values is now 576) i.e. the relative yield change is only 
7%. However, the change in the hyperfine coupling is 
about 300%, clearly exemplifying the disparity of clas- 
sical expectations with observed data. We will now fo- 
cus on more recent data regarding the effect deuteration 
has on the magnetic sensitivity [sol . [sij of the the pyrene 
(Py)-dimethylaniline (DMA) RIP recombination dynam- 
ics. What has been observed is that the low-field-effect 
is most pronounced when only one of the radicals ex- 
hibits strong hyperfine couplings, with the second radi- 
cal being deuterated and hence having substantially re- 
duced hyperfine couplings (the magnetic moment of deu- 
terium is about one third of the proton's). In particular, 
in the case of protonated pyrene RIPs Py-h^g/DMA-h^^j^ 
and Py-hj'p/DMA-d;^';^, there was almost no magnetic- 
sensitivity of the exciplex fluorescence. On the contrary, 
the recombination dynamics of deuterated pyrene RIPs 
Py-dJ"o/DMA-h:J';^ and Py-dj'o/DMA-dfi showed a clear 
low-field sensitivity. In Figure 6 we show the calculated 
behavior of the magnetic sensitivity. For each of the four 
RIPs we calculate the eigenvalue spectrum and count the 
decay rates A that are smaller than lo. This spectrum is 
shown for one of the RIPs in Figure 6A. Half the percent- 
age of those decay rates is a good estimate 26] for the 
decrease in the singlet recombination yield. It is seen that 
the relative yield change with magnetic field displayed in 
Figure 6B is in perfect agreement with observations. The 
absolute values of the yield drop are not accurately pre- 
dicted due to the simplicity of the two-spin model used 
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FIG. 6: Magnetic sensitivity of the reaction yield of the Py- 
DMA radical-ion pair. What is calculated is the drop in 
the singlet-recombination yield. The magnetic Hamiltonian 
contains two nuclear spins, with Ji and I2 being either 1/2 
or 1, dependent on which of the four RIPs is considered. 
The hyperfine couplings are taken from [3(1 |. The singlet 
and triplet recombination rates used were ks — 50 /is~^ and 
kr = 0.2fcs, respectively. Similar results where obtained for 
ks = 10,20 ^s-^ (A) The decay rates of Py-h^o/DMA-h+ 
as a function of magnetic field. Similar plots are obtained 
for the other radical-ion pairs. The magnetic sensitivity is 
found by counting the decay rates that are smaller than uj. 
(B) Estimated singlet-yield drop from the count of the slow 
decay rates for all four Py-DMA RIPs. It is seen that the 
protonated pyrene pairs show no magnetic response. 



(they also depend on the particular value oi ks). 

In summary, we have analyzed the physical interpreta- 
tion of low magnetic field effects that is bestowed upon 
the radical-ion pair recombination reactions by their 
full quantum-mechanical description. The explanatory 
power of this theory has been applied in understanding 
experimental observations on the magnetic sensitivity of 
RIP reactions that were conflicting classical reaction the- 
ories. It has been shown that the eigenvalue spectrum of 
the density matrix equation describing the time evolu- 
tion of the rip's spin state accounts for the observations 
regarding RIP deuteration. In particular, the slow decay 
rates imposed by the quantum Zeno effect on the eigen- 
value spectrum, and their dependence on the problem's 
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physical parameters, are shown to contain rich dynamical effect of low magnetic fields on photosynthetic reactions 
information. Finally, it is rather crucial to address the under this new perspective. 
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